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The paper contains some new results and a review of recent achievements, concerning the multisup- 
■ port solutions to matrix models. In the leading order of the 't Hooft expansion for matrix integral, 

these solutions are described by quasiclassical or generalized Whitham hierarchies and are directly 
related to the superpotentials of four-dimensional Af = 1 SUSY gauge theories. We study the 
derivatives of tau- functions for these solutions, associated with the families of Ricmann surfaces 
(with possible double points), and relations for these derivatives imposed by complex geometry, 
including the WDVV equations. We also find the free energy in subleading order of the 't Hooft 
expansion and prove that it satisfies certain determinant relations. 



O 

in 



Q_i| Recent interest to matrix models and especially to their so-called multisupport (multicut) solutions 

was inspired by the studies in J\f = 1 SUSY gauge theories due to Cachazo, Intrilligator and Vafa ^P, 
[2] and by the proposal of Dijkgraaf and Vafa [Sj to calculate the low energy superpotentials, using 
the partition function of multicut solutions. The solutions themselves are well-known already for a 
long time (see, e.g., [HE]) with a new vim due to the paper by Bonnet, David and Eynard [§]. 

H ' , 

The Dijkgraaf- Vafa proposal was to consider the nonperturbative superpotentials of N = 1 SUSY 

gauge theories in four dimensions (possibly coming as the softly broken J\f = 2 Seiberg-Witten (SW) 

theories (Z1IH1) arising from the partition functions of the one-matrix model (1MM) in the leading order 

in 1/N, N being the matrix size. The leading order (of the 't Hooft 1/iV-expansion) of the matrix 

model is described by the quasiclassical tau-function of the so-called universal Whitham hierarchy 9 

(see also the details about one-matrix and two-matrix cases see in ^5] and [TH]b 

One would expect the existence of the relation between the partition function in the planar limit and 

the quasiclassical, or Whitham hierarchy already because matrix integrals are tau-functions of the 

hierarchies of integrable equations of the KP/Toda type [T7] . For the planar single-cut solutions the 

matrix model, partition functions become tau-functions of the dispersionless Toda hierarchy, one of 

the simplest example of the Whitham hierarchy. 

One may also consider more general solutions to matrix models, identifying them with generic 
solutions to the loop (Schwinger-Dyson, or Virasoro) equations JH]) to be the Ward identities satisfied 
by matrix integrals These generalized solutions to the loop equations can be still treated as 
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matrix model partition functions (e.g., within a D-module ideology |2L)l 1211 l2~2*] l: however, they do not 
necessarily admit any matrix integral representation. 

In what follows, we aim to discuss an interesting class of multi-cut, or multi-support, solutions to 
the loop equations that have multi- matrix integral representation I23| I20| I22j . These solutions are 
associated with families of Riemann surfaces and form a sort of a basis in the space of all solutions to 
the loop equations |201 122j (like the finite-gap solutions form a basis in the space of all solutions to an 
integrable hierarchy). They can be distinguished by their "isomonodromic" properties — switching on 
higher matrix model couplings, or 1/iV-corrections does not change the family of Riemann surfaces, 
but just reparameterizes the moduli as functions of these couplings. This property is directly related 
to that the partition functions of these solutions are quasiclassical tau-functions (also often called as 
prepotentials of the corresponding Seiberg-Witten-like systems). 

In sect.^ we describe the general properties of multi-cut solutions of matrix models. In particular, 
we prove that the free energy of the 1MM in the planar (large N) limit coincides with the prepotential 
of some Seiberg-Witten-like theory. This free energy is the logarithm of a quasiclassical tau-function. 
The corresponding quasiclassical hierarchy is explicitly constructed. 

The Whitham hierarchy is basically formulated in terms of Abelian differentials on a family of Rie- 
mann surfaces [S]. This implies the main quantities in matrix models are to be expressed in geometric 
terms and allow calculating derivatives of the matrix model free energy. Indeed, we demonstrate in 
sect. 121 that the second derivatives of the logarithm of the matrix model partition function can be 
expressed through the so-called Bergmann bi-differential. 

In sect. |21 we turn to the third derivatives of partition function and to the Witten-Dijkgraaf- 
Verlinde-Verlinde (WDVV) equations [241 1251 l2*S] . which are differential equations, involving the third 
derivatives of tau-function with respect to Whitham times. These equations are usually considered an 
evidence for existing an underlying topological string theory. In sect. 3, we prove that the quasiclassical 
tau-function of the multi-support solutions to matrix models satisfies the WDVV equations in the case 
of general 1MM solution, i.e., in the case of arbitrary number of nonzero times, which include now, 
besides the times of the potential, the occupation numbers (filling fractions) indicating the portions 
of eigenvalues of the corresponding model that dwell on the related intervals of eigenvalue supports. 
Although being, at first glance, quantities of very different nature in comparison with the original 
times, they can be nicely combined into a unified set of the "small phase space" of the model 1 . This 
completes an interpretation of the results of pj |2j in terms of quasiclassical hierarchies. 

The WDVV equations are a simple consequence of the residue formula and associativity of some 
algebra (e.g., of the holomorphic differentials on the Riemann surface) |26[ I27| I28j (see also [29 J); 
moreover, the associativity of algebra can be replaced by a simple counting argument [20]: the number 
of critical points in the residue formula should be equal to the dimension of the "small phase space". 
We present the proof of the residue formula and extra conditions following |31j . Note that, while 
the residue formula is always present in the theories of such type [2H]) 2 the associativity (naive 
counting) is often violated 28 ( 30 ) (maybe the most notorious case is the SW system associated 
with the elliptic Calogero-Moser model, where this phenomenon was first found in |28|). 

There are strong indications that the correspondence between matrix models and SUSY gauge 
theories goes beyond just the large- N limit of matrix models. Say, the relation between gauge theories 
and matrix model were further verified at the nonplanar, genus-one level for the solution with two 
cuts and a cubic matrix model potential [221- In sect. [IJ we find the multicut solution to 1MM in 
the subleading order: the torus approximation in terms of a dual string theory. We calculate the 
corresponding free energy and, in particular, prove it to have a determinant form, related to the 

1 Note that, from the point of view of N — 1 SUSY theory, couplings in the matrix model potential must be identified 
with couplings in the tree superpotential, while the occupation numbers play the role of moduli of the loop equation 
solutions and must be associated with the vacuum expectation values of the gluino condensates. 

2 The recently proposed in |32l residue formulas in planar case just follow from the standard residue formula. However, 
the residue formulas of |32| involving non-planar corrections look new and rather instructive. 
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topological B-model on the local geometry II and conjectured in .'5 (the authors of 2.'} have checked 
it for the several first terms of expansion in the two-cut case). Then, we apply formulas for the third 
derivatives of the planar free energy and for the genus-one free energy obtained in the paper, and 
conjecture a diagrammatic interpretation of these formulas to be extended to higher genera. 



1 Matrix models and generalized Whitham hierarchies 

1.1 Matrix integrals and resolvents 

Consider the 1MM integral 3 

/ DXe~^ v ^ =e r , (1) 

J NxN 

where V(X) = J2 n >i^nX n , Ti = is a formal expansion parameter, the integration goes over the 
NxN matrices, DX oc Yiij dXij, and for generic potential one should consider the holomorphic version 
of (^P) implying contour integration in complex plane in each variable. The topological expansion of 
the Feynman diagrams series is then equivalent to the expansion in even powers of U for 

oo 

T = to,h,h, • • •) = E h 2h ~ 2 Fh, (2) 

h=0 

Customarily to = tiN is the scaled number of eigenvalues. We assume the potential V{p) to be a 
polynomial of the fixed degree m + 1, with the fixed constant "highest" time t m +i- 

The averages, corresponding to the partition function JQ) are defined as usual: 

(f(X)) = \\ DXf(X) exp f-hvV(X)) (3) 
6 J NxN \ n J 



and it is convenient to use their generating functionals: the one-point resolvent 

trX* 



mx) = n± { -^ (4) 



k=0 

as well as the s-point resolvents (s > 2) 



fe, =1 • • • As \ A l A s -a. I conn 



k\ ,...,k 3 =l 

where the subscript "conn" pertains to the connected part. 

These resolvents are obtained from the free energy T through the action 



W(Ai,...,A.) = ti 



2 d d or 



dV{\ s ) dV(\s-i) 9V{\ X 
d d d 



W{\i), (6) 



dV(X s ) dV(X s -i) dV(X 2 ) 

3 There is a more consistent point of view on matrix integrals and loop equations below that implies introducing 
some background (polynomial) potential Vo(x) into the exponential Q, i.e., shifting first several tk — > Tk + tk and then 
looking at this matrix integral as a formal series in t^s (see, e.g., |33l I2UI for a review). Then, the loop equations 
are iteratively solved. Such a framework is more effective when looking at the whole manifold of solutions to the loop 
equations. However, this is at the price of rigid fixing the number of cuts in the multicut solution from the very beginning. 
Since we will freely change the number of cuts (which is a smooth procedure specifically in the multicut solutions) and 
do not care of other than multicut solutions to the loop equations, we follow a less pure scheme (sec, e.g., 4 6 34 and 
references therein) that basically would imply some re-summation of infinite series, etc. It allows us, however, instead of 
dealing with infinite series to deal with objects determined on Riemann surfaces. 



3 



of the loop insertion operator 4 

dV(X) dtj v ; 

Therefore, if one knows exactly the one-point resolvent for arbitrary potential, all multi-point resolvents 
can be calculated by induction. In the above normalization, the genus expansion has the form 

oo 

W(\ l ,...,\ s ) = Y / h 2h W h (\i,...,\s), s>l, (8) 

h=0 

which is analogous to genus expansion 

The first in the chain of the loop equations ^H] of the 1MM is ^5] 

f d^Y!i^l W (\) = KW(x) = W(x) 2 + h 2 W(x,x) (9) 
Jct, 27ri x — X 

where the linear integral operator K, 

projects onto the negative powers 5 of A. Hereafter, Cx> is a contour encircling all singular points of 
W(X), but not the point A = p. Using Eq. (JHJ), one can express the second term in the r.h.s. of loop 
equation © through W(p), and Eq. © becomes an equation for the one-point resolvent Q. 

Substituting the genus expansion (|SJ) in Eq. @, one finds that Wh(X) for h > 1 satisfy the equation 

(K- 2Wo(A)) W h (X) = ]T W h r(\)W h - h r(\) + ^— rf W (A), (11) 



h'=l 



In Eq. 1)11(1 . Wh{\) is expressed through only the Wh(X) for which /ij < /i. This fact allows one to 
develop the iterative procedure. 

To this end, one needs to use the asymptotics condition (which follows from the definition of the 
matrix integral) 

W h (X)\ x ^oo = jS h ,o + 0(1/A 2 ), (12) 

and manifestly solve © for genus zero. Then, one could iteratively find Wh(X) and then restore the 
corresponding contributions into the free energy by integration, since 

Wh{x) = Wixf^ h - 1 ' (13) 

1.2 Solution in genus zero 

In genus zero, the loop equation ® reduces to 

dX V'(X) 



Ct, 27TZ X — X 



W (X) = [V'{x)f{x)\ _ = (W (x)Y (14) 



4 This operator contains all partial derivatives w.r.t. the variables t^'s. However, below we introduce additional 
variables 5; and, therefore, we use the partial derivative notation here. 

5 In order to prove it, one suffices to deform the integration contour to infinity to obtain 

^^A /(A) = V'^fW ~ [V'(x)f(x)] + = [V'(x)f(x)}_ 

C-D 
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In order to solve this equation for the planar one-point resolvent Wo(p), one suffices to note that 

[V(X)W (X)]_ = V'(X)W Q (X) - [V'(X)W (X)] + (15) 

and, due to (J12j) . the last term in the r.h.s. is a polynomial of degree m — 1, m being the degree of 
V'(X), 

P m -iW = -[V'(X)W (X)] + = -f c ^-YM Wo (x) (16) 
It can be also rewritten as action of the linear operator ry(A) acting to the free energy, 

P m -i(A) = -ry(A)^ , MA) = E( fe + 1 + 2)^+ 2 A fc ^- (17) 



Then, the solution to (|14[) is 



1 



Wo(A) = -V'{X) - -yJV>(\)* + 4P m _ 1 (A), (18) 

where the minus sign is chosen in order to fulfill the asymptotics (|12(l . For the polynomial potential 
of power m + 1, the resolvent VFo(A) is a function on complex plane with m cuts, or on a hyperelliptic 
curve 

y 2 = V'(Xf + 4P m _x(A) (19) 
of genus g = m — 1, conveniently presented introducing new variable y by 

W (X) = ±{V'(X)-y), (20) 

For generic potential V^(A) with m — > oo, curve (|19|) may have an infinite genus, but we can still 
consider solutions with only finite number n of cuts and separate the smooth part of curve (|19|) 
introducing 

y = M(X)y, and "reduced" Riemann surface y 2 = ]^[ _ 1 (A — fJ. a ) (21) 

with all iJL a distinct. In what follows, we still assume M(A) to be a polynomial of degree m — n, keeping 
in mind that n is always finite and fixed, while m > n can be chosen arbitrarily large. By convention, 
we set j/Ia-^oo ~ A n , and M(A) is then 6 

r A V 1 ! \ m—n 

w-L&v^mr)*"^-" (23> 

Note that the values of M(A) at branching points, M m _ n (ju a ) coincide with the first moments Ma^ 
of the general matrix model potential (see |35U34l |2"3] ) . while the higher moments are just derivatives 
at these points. Indeed, by their definition, these moments are 

M^=i — ( n f , r , p>l (24) 
lev 2iti (x - ix a )Py(x) P ~ V ' 

Then, one deforms the integration contour to infinity assuming V'{x) is an entire function and, as 
in (|23|) . make use of the asymptotic conditions (|22j) . to replace V'(x) in the numerator of (|24l) by 
y(x) = M m - n {x)y{x) subsequently evaluating the integral by taking the residue at the point x = fi a 
and obtaining 

1 / d p ~ l 



, e.g., Ml') = M{ix a ). (25) 
A=/i Q 



3 Since, due to 1191 . at large A 



\/'(A)=y(A) + ^i + 0(l/A 2 ) (22) 
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Inserting the solution (|21j) . 1)23(1 into l|2U|) and deforming the contour, one obtains the planar 
one-point resolvent with an n-cut structure, 



WoW =y ±mm (26) 

2 Jc-d 27tz A - x y(x) 

The contour Cx> of integration here encircles the finite number n of disjoint intervals 

n 

V = [J\P2i-U V>2i]> A*l < A*2 < • • • < M2n- (27) 
i=l 

Let us now discuss how many free parameters we have in our solution. If one does not keep genus 
of the curve fixed, it is given for a generic potential V(X) by m + 1 times t&, coefficients of V'(A), and 
m — 1 coefficients p& of the polynomial P m _i(A), its leading coefficient being related to to- Indeed, 
using (|T2|) and (fl8|) . one obtains 

««A)lv~ = "%^f + ■ ■ ■ - " (^g^ + °< VA 2 ) = | + °(1/A") (28) 

Therefore, totally one has 2m + 1 parameters (including to)- in fact, we usually fix the leading 
coefficient of the potential V(A) to be 1/m (see s.3), i.e. M = 1 in l|23j) . which leaves us with 2m 
parameters. 

If, however, now one fixes the curve of genus g = n— 1, (|21|) . this imposes m—n conditions of double 
points (coinciding branching cuts) 7 . Therefore, one then has 2m — (m — n) = m + n parameters. Of 
these, still m + 1 parameters are tk, which are variables in the loop equations, while n — 1 parameters 
give the arbitrariness of solutions to these loop equations (possible different choices of the polynomial 
Pm-l)- 

One can arbitrary choose coordinates on this n — 1-dimensional space of parameters. It turns out, 
however, that there is a distinguished set of n independent variables that parameterize solutions to 
the loop equations El , 

Si=f ^V={ ^M(X)y : (30) 



\Ai 4vrz J a, 47rz 

where Ai, i = 1, . . . , n — 1 is the basis of A-cycles on the reduced hyperelliptic Riemann surface (|21[) 
(we may conveniently choose them to be the first n — 1 cuts) see Fig.l. 

Besides canonically conjugated A- and -B-cycles, we also use the linear combination of 5-cycles: 
Bi = Bi — Bi + i, B n -\ = B n -\. Therefore, -B-cycles encircle the nearest ends of two neighbor cuts, 
while all -B-cycles goes from a given right end of the cut to the last, nth cut. For the sake of definiteness, 
we order all points \x a in accordance with their index so that fi a is to the right of fxp if a > j3. 

1.3 Matrix eigenvalue picture: a detour 

The variables Si can be formally determined as eigenvalues of a differential operator in t&'s |22j . How- 
ever, they can be also more "physically" interpreted in the quasi-classical picture of matrix eigenvalues 
(Coulomb gas) when their number (and, therefore, size of the matrix) goes to infinity. We come now 
to this interpretation. 

7 One can also give these conditions via different integral conditions. Say, one can get from 1121 and 12611 . the 
asymptotic conditions 

1 / d\ \ k V'(X) , „ 
t Sk,n = -<p — k = 0,...,n. 29 

These conditions are identically satisfied for m — n, since V'(X) ~ J/(A) + 0(1/A) = j/(A) +0(1/A), while for m> n they 
impose constraints. 



6 



Si S2 S3 Sg t() — J2Si 

A, B l A 2 ^ ^ ^ Bg-Bg ^ 



Figure 1: Structure of cuts and contours for the reduced Riemann surface. 

To this end, let us first introduce the averaged eigenvalue distribution 

t N 1 
^ E W A " *0> = T7- lim (w(X - ie) - W(X + iej) (31) 



N ^ x y " 2m €^0 
where Xj's are eigenvalues of the matrix X. In the planar limit, this quantity becomes 

Po (\) = -Llim(Wo(\-ie)-Wo(\ + ie)) = ^-lmy(X) (32) 

and satisfies the equation 8 

/ = \v'{x), Vp G V (33) 

Jv x — X 2 

This averaged eigenvalue distribution becomes the distribution of eigenvalues in the limit when their 
number goes to infinity. For the illustrative purposes, let us do the matrix integral, ((JJ performing 
the matrix X in the form of U • diag(xj) • U~ l with a unitary matrix U and then first making the 
integration over the "angular" variables U. Then, one comes to |36| 

/n^iii^-^) 2e "^ Eiy(a!i) = 

i i>j 

= I JJ (icfe -^(/nA)£»(A)-/ e (A) <f (A')iog|A-A'|d^O = J Yldxte^ff (34) 

% i 

where we introduced the eigenvalue distribution 

<?M = | ( 35 ) 

i 

Now, in the limit of large N, one can use the saddle point approximation to obtain the equation for 
g(X). However, one also need to take into account the constraint 



g(X)dX = t (36) 

by adding to S e //, the Lagrange multiplier term IIq(/ g—to). Note also that g(X) is a non-negative 
density. This finally leads to the saddle point equation 

2 J g(X) log \x - X\dX = V(x) + U , \/x G support of g (37) 

The derivative of this equation w.r.t. p coincides with (|33j) . 



'indeed, by definition 



W ^-dX = 



x — A 

Now, using [12( )i and the definition 1311 1 and pulling out the contour from infinity, one easily comes to this equation. 
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The variable to plays role of the (normalized) total number of eigenvalues, 



If 1 

to = 4tH Jc = ~2 VeSoc ( yd ^ ( 38 ) 

and the support of g is T> that consists of N segments T>i. Now, following [HJ one may fix the 
(occupation) numbers of eigenvalues in each of the segments, Si (|3U|) . i = 1, . . . ,n — 1. We assume the 
occupation number for the last, nth cut to be to — J27=i Si = S n . 9 (Obviously, no new parameters 
Si arise in the one-cut case.) We formally attain this by introducing the corresponding chemical 
potentials (Lagrange multipliers) Hi, i = 1, . . . , n — 1, in the variational problem for the free energy, 
which therefore becomes in the planar limit 

S eff [Q;Si,to,t k ) = -J V{\) Q {\)d\ + J J Q(X)\og\X-X'\ e {\')dXdX' + 

n-l 



-n Q (X)dx - 1 ) - E n * (J D ^ x ) dX ~ s i) ■ ( 39 ) 

while the saddle point equation becomes 

2 J q{X) log \x - X\dX = V(x) +Ui + n , Vx E V { (40) 

and its derivative still coincides with (|33|) . 

Therefore, with generic values of the constants IT, £ C (A) gives the general solution to (JHH]) (or the 
planar limit of the loop equation) : these constants describe the freedom one has when solving the loop 
equation 10 . However, in the matrix model integral (where there are no any chemical potentials) one 
would further vary To w.r.t. IT to find the "true" minimum of the eigenvalue configuration, 

This is a set of equation that fixes concrete values of Si and IT in the matrix integral. 
Let us now calculate the derivative of To w.r.t. Si. From (|39|) . one has 

dSeff 



dSi 



Q=Pu 



V 



dA ^S~ ( y(A) - 2 f v dX ' lo ^ X ~ A >(A')) (42) 



The expression in the brackets on the r.h.s. of (|42[) is almost a variation of (|«S9j) w.r.t. the eigenvalue 
density, which is 



_ 5S eff 



8p(X) 



Q=P0 



V(X) -2 [ dX' log(A - A>o(A') + V-i + n 
Jv 



for A G A C V. (43) 



It is therefore a step function, h{X) which is constant equal to Q = — IIo — H on each cut Ai. One 
then has 



dTo _ dS eff [g] 



dSi OS, 



d 

OS 



EW = -6 + 6i = ^ (44) 



9 It is sometimes convenient to consider instead of to as a canonical variable. However, in all instants we use S n , 
we specially indicate it for not confusing S n with the "genuine" filling fraction variables Si, i = 1, . . . , n — 1. 

10 Note that, instead of fixing the occupation numbers, one could use other ways to fix a solution to the loop equations, 



sec 
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In particular, 

f - n « <«> 

In [I] it was proved that the difference of values of ITj on two neighbour cuts is equal to 11 

O+i - Ci = 2 / W (A)dA (46) 

i.e. 

IT = (Ci+i - Ci) + (Ci+2 - + • • • + (Cn-i - Cn-2) + (Cn - Cn-i) = / yd\ = I ydX (47) 

•/B i UB i+ iU...UB 9 iB, 

Note that one can calculate the planar limit free energy that can be obtained via substituting the 
saddle point solution g into (f3*9l is 

n— 1 

= s e // [&] = -\ I v(\) 6c (\)d\ + -n i + n * 5 * ( 48 ) 

In the paper, we choose the solution to the loop equation with fixed occupation numbers, (|3U[). 
Note that fixing the chemical potentials (|44| ) -(|47| ) instead, (|41|l. corresponds just to interchanging 
A- and -B-cycles on the Riemann surface (|19|l . However, is n °t modular invariant. Under the 
change of homology basis, J-q transforms in accordance with the duality transformations [HZj (which 
is a particular case of behaviour of J-q under the general transformations, 38 ). The higher-genus 
corrections become also basis-dependent: choosing Si or IT as independent variables, one obtains 
different expressions, say, for the genus-one free energy, see s.4.3. 

In the next two subsections we are going to demonstrate that the planar loop equation solution 
with fixed occupation numbers corresponds to a Seiberg-Witten-Whitham system. 



1.4 Seiberg— Witten-Whitham theory 

Seiberg-Witten system. We call the SW system |7f 12 the following set of data: 

• a family A4 of Riemann surfaces (complex curves) C so that the dimension of moduli space of 
this family coincides with the genus 13 ; 

• a meromorphic differential dS whose variations w.r.t. moduli of curves are holomorphic (this 
implies existence of a connection on moduli space, so that this statement has a strict sense, see 

e.g. nam]). 

These data allow to define a SW prepotential [21 related to an integrable system EH IS El EQ| • 
First, one introduces the variables (whose number coincides with the genus of C) 

Si = — I dS (49) 

where Ai are A-cycles on C. As soon as |rgr is holomorphic, from the definition l|49j) of Si and the 
obvious relation dSj/dSi = 5ji, one finds that 

1 r ddS r 

6 jh (50) 



4m 7a, dSi 



lx The simplest way to prove it is to define function h(\) outside the cuts: h(X) — V(A) — 2 J v dA'log(A — A')po(A') 
and note that ^'(A)^^ = 2W (X). 

12 Various properties of such systems can be found in |12l 1131 HI] . 

13 More generally, this can be extended to the curves with auxiliary involution, or certain directions in moduli space 
can be "frozen" in a different way, certain examples can be found in |28|. 
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i.e., 

-as; = d ^ (51) 

where dcoi are the canonically normalized holomorphic 1-differentials, 

±-J a ^ 3 = 5 13 (52) 

Introducing £>-cycles conjugated to ^4-cycles: Ai o Bj = 5(j, where o means intersection form, one 
obtains that 

^sJ Bj dS= i^ =Tii (53) 

is the period matrix of C and is therefore symmetric 14 . Hence, there exists a locally defined function 
F such that 

wri iS - <55) 



and this function is called a prepotential. 



Generalized Whitham system. We now extend the Seiberg-Witten system to the quasiclassical 
or generalized Whitham system [:9; by introducing extra parameters or times t k into the game, we do 
it mostly following |10 | lllj 15 . In order to construct a Whitham system, one needs to add to the SW 
data a set of jets of local coordinates in the vicinity of punctures on C. In particular case of (|19|) . 
1)21(1 and these points (the singularities of ydX coincide with two A-infinities, where we choose the local 
parameter £ = j. We then introduce a set of meromorphic differentials dQ. k with the poles only at 
these punctures (as the hyperelliptic curve (|19|) . (|21|) is invariant w.r.t. the involution y — > —y, from 
now on we just work with either of the two infinities, see |10| lllj) with the behavior 

dn k = ±| (V fc-1 + O(l)) d£, for £ 0, k > (56) 

where signs are different for the different infinities. We also introduce the bipole or 3rd-kind Abelian 
differential, with two simple poles, located at two infinities: 

reSoodfio = — resoo_(ir2o = — 1 (57) 

Then, the generalized Whitham system is generated by a set of equations on these differentials and 
on the holomorphic differentials dujc 

dt k = dt p ' dSi = ~dh' ~dS~ = ~dS~ ^ ' 

where the partial derivatives are supposed to be taken at constant hyperelliptic co-ordinate A. 
These equations imply that there exists a differential dS such that (see also ()51(l 

ddS , ddS ln . . 

— = dia l , —- = 2dn k 59 
dSi dt k 



'This follows from the Riemann bilinear relations for canonical holomorphic differentials 15111 
= / duii A duij — 

k 

= 47rz(T ij - 7; 7i ) (54) 



15 One can introduce these extra parameters in the context of supersymmetric SW gauge theories, even without reference 
to an integrable system a priori, see |42j and references therein. 
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Note, however, that the meromorphic differentials (|56|) . ()57() are defined up to linear combinations of 
the holomorphic differentials. Since we consider Si and i/% to be independent variables, this ambiguity 
is removed merely by imposing the condition |1UI lllj 

OS 1 P 

—i = — J> dn k = ovi,k (60) 

dt k 2m J a, 

Now we can invariantly introduce variables t k via the relations 16 

t k = --^resoo (\~ k dS^J , k = 1, . . . , m, t = res^dS (cf. ©) (61) 

Then, one defines the prepotential that depends on both Si and t k via the old relation (J55|) and the 
similar relations 

OF 1 1 

— = ~res\ =00 (\ k dS) = ~v k , k=l,...,m. (62) 
OF r°°+ 

arrL dS <63) 

The latter integral, which is naively divergent, is still to be supplemented with some proper regular- 
ization, we discuss this in detail in the next subsection. 

In fact, one still needs to prove such a prepotential exists jlUI II lj . by checking the second derivatives 
are symmetric. This can be verified using the Riemann bilinear relations, 



^ = 2reSoo (Q k duJi) = I 
obi Jdi 



Q, k duJi 



2 Y(I ^tdui- I nzduj l )-2y ( f nUui- ( VLZdwi 

r-f KJBi Jb 1 ) r-r \ja 1 ja 1 



i=i ud i ' i=i 



2^ ( <z> cZrifc f diOi — f dflk f duji 1=2® dVL k = — — . (64) 
KJa, JBi JBi J A t J Jb, ot k 



Here 0T, g is the cut reduced Riemann surface (|21j) (see fig. |2J) , f2j^ are values of Q/. = J dfl k on two 
sides of the corresponding cycle (they are related by the integral over the dual cycle), and we have 
used (J72J). 

Similarly, for the derivatives w.r.t. the times tk, one has analogously to (|64|) 



1 / dvk_ _ dvp_ 

2 I dt p dt k 



resoo ((Q k ) + d£l p - dfl k (fl p ). 
resoo (QkdQ, p ) = <z> fifcdfL 

JdT, a 



l=1 \JAi JB t JB t J Ax 



2_. ( f dVL k f dQ p — 6 dOfc f dtt p I = 0. (65) 



and 



= / duji A dQ,Q = ( f dtOi f dClo — f dQo f duj>, 
JY; g ^ \JA t JBi J A, JBt 

+res 00 (du>i) / dQo — ies 00 (dQ,o) / dw{ = 

Joo— Joo— 



+ 



noo+ 

dQ,Q — I du>i. (66) 

Bi J oo— 



16 In what follows, we call the oo-point the point oo + , or A = oo on the "upper" (physical) sheet of hyperelliptic 
Riemann surface (I2H corresponding to the positive sign of the square root. 
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Figure 2: Cut Riemann surface from Fig. El The integral over the boundary can be divided into several 
pieces (see formula (|118)l ). In the process of computation we use the fact that the boundary values of 
Abelian integrals lo^ on two copies of the cut differ by period integral of the corresponding differential 
duj over the dual cycle. We add the two infinity points and the additional, nth, cut dividing the 
surface in two sheets. We present the logarithmic cut between these two points (about which we draw 
the standard integration contour C£). Integrals over small circles around the points u a are relevant 
only when calculating the third-order derivatives w.r.t. the canonical variables t/. 



1.5 Free energy as prepotential of SWW system 

We now associate a Seiberg-Witten-Whitham (SWW) system with the planar limit of the matrix- 
model free energy. 



Matrix integral as a SW system. The family Ai in this case is the family of h = n — 1 reduced 
Riemann surfaces described by (|19|) or (|21j) . In different words, there is no information in Ai about 
the additional polynomial M m _ n (A), which is present, however, for (|21|) in the differential dS. The 
role SW differential is played by 

dS = ydX (67) 

Consider its variation w.r.t. Si, the variation over moduli Ai does not change the genus of the reduced 
Riemann surface as well as the highest degree of the additional polynomial M m _ n (A). Moreover, 
considering the times of the potential V^(A) to be independent on the parameters Si, we assume 
5V'/SSi = 0. Below, by 5 and 5s we denote the respective general variation and variation specifically 
w.r.t. the moduli parameters Si. 

Using (j!9j) . (|21|) . one obtains for the general variation 5dS: 17 



5dS = 5 (M m _ n (A)y(A)) dX 

iM (\)5n^(X) 

-dX, (68) 



#2n(A)<5M m „ n (A) + \M m - n (X)5g 2 n{X) 



m 

where the polynomial expression in the numerator is of maximum degree m + n — 1 (since the highest 
term of M m _ ra is fixed). On the other hand, under 5s which does not alter the potential, we obtain 
from (UU), (JUJ) that 

r jq 1 5 S P m ~l{X) 

5 s dS = --— m - m Q*- (69) 

2 M m _ n (A)y(A) 



17 Note that the variation S differs nevertheless from loop insertion (J7J because the former does not change, by definition, 
the degree of the polynomial Af m _ n (A). 
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Because this variation is just a particular case of l|68|l. we obtain that zeros of M m _ n (A) in the 
denominator of ()69|) must be exactly cancelled by zeros of the polynomial 5sP m ~i (A) in the numerator, 
so the maximum degree of the polynomial in the numerator is n — 2 (because, again, the highest-order 
term of P m _i(A) is fixed by asymptotic condition (f2*2*|) and is not altered by variations 5s)- We then 
come to the crucial observation that the variation SgdS is holomorphic on the curve (|21j). as it should 
be for the SW differential. 

The canonical 1-differentials on the reduced Riemann surface y(X) have the form 

ddS Hj{\)d\ 

-twr dUi = ~wr (70) 

and -ffj(A) are polynomials of degrees at most n — 2. The normalization condition (|52|) unambiguously 
fixes the form of the polynomials Hi(X). 

Comparing now (|47|) and (|55jl. one could indeed identify the planar limit of the 1MM free 
energy with an SW prepotential F. However, in order to specify this equivalence further, one needs 
to work out the t-dependence of the free energy, i.e. consider the generalized Whitham system. 



Matrix integral as a Whitham system. To this end, let us check that differential (|fi7j) on curve 
(JUJ with the relation for moduli lfT9"|). (J2TJ) does satisfy 

Indeed, we have proved the first set of relations (|59j) in the previous paragraph. Now let us consider 
variations of the potential, i.e., variations w.r.t. Whitham times t k . Then, we obtain instead of ()69j) 

2 M m _„(A)ti(A) ' 

while (|68[) still holds. Repeating the argument of the previous paragraph, we conclude that the 
zeroes of M n ^ k (X) cancel from the denominator and, therefore, the variation may have pole only at 
A = oo, or rj = 0, i.e., at the puncture. In order to find this pole, we use ()68|) . which implies that 
dS = M m _ n (A)y(A)dA — > (V^ n (X) + 0(j))dX and, therefore, the variation of dS at large A is completely 
determined by the variation of V^X). Parameterizing V(X) = J2T=i^kX k , we obtain up to a 
linear combination of holomorphic differentials. One may fix the normalization of d£lk that are also 
defined up to a linear combination of holomorphic differentials in order to make Eq. (|59|) exact. We 
already discussed (see lfiU|) that this normalization and, therefore, unambiguous way the variables Si 
depend on the coefficients of P m -\ is fixed by the condition [TU1 ITT] 

dS- 1 f dS- 

afr^L dnt=0V '- k - or 8#r° i = i --»- 1 < 72 > 

The derivatives of dS w.r.t. the times are 

njn _ ddS V'(X)kX k - 1 dX 1 ^ 2 dP, X j dX 
2dU k = — = + 



dt k y 2 dt k y 

H n+ k-i{X)dX 



(73) 



and the normalization conditions (|72j) together with the asymptotic expansion 

2dn k {X)\ X ->oo = kX k - 1 dX + 0(X- 2 )dX = 

oo 

= kX k ~ l dX + c fcm A -1 - m dA. (74) 

m=l 

fixes uniquely the coefficients of the corresponding polynomials H n+k _\ of degrees n + k — 1. 
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Now we again find that the prepotential F defined in Q62|) coincides with Tq. To this end, we 
apply the formula similar to 1)42 jl with the only difference that the potential V(A) itself is changed. 
We then obtain (see flHJ)-(|IIJ)) 



3^0 



dt k 




(75) 



which by virtue of (|72"|) gives (|62|). 

On to-dependence of prepotential. Thus, we have proved the derivatives of the SWW prepoten- 
tial and of the matrix model free energy w.r.t. Si's and tk's coincide. However, there is a subtlety 
of exact definition of integral Q63|) for dF/dtQ. This quantity is to be compared with the derivative 
dFo/dto, ({I!)) equal to the integral § v log |A — y\dS — V(y), where the reference point y is to be chosen 
on the last, nth, cut, while the expression itself does not depend on the actual local position of the 
reference point. It is convenient to choose it to be fi2n = b n — the rightmost point of the cut. We can 
then invert the contour integration over the support D to the integral along the contour that runs first 
along the upper side of the logarithmic cut from b n to a regularization point A, then over the circle 
C\ of large radius |A| and then back over the lower side of the logarithmic cut in the complex plane. 
In order to close the contour on the hyperelliptic Riemann surface under consideration, we must add 
the integration over the corresponding contour on the second sheet of the surface as shown in Fig. 01 
we let Cl denote the completed integration contour, and it is easy to see that such an additional 
integration just double the value of the integral. 

It is easy to see that all the singularities appearing at the upper integration limit (i.e., at the 
point A) are exactly cancelled by the contribution coming when integrating the expression dS log(A — 
b n ) along the circle C\; in fact, the latter can be easily done, the result is — 2-7ri(5(A) — {S(b n )} + ), 
where the function S{\) is the (formal) primitive of dS (which includes the logarithmic term), and the 
symbol {•}+ denotes the projection to the strictly polynomial part of the expression in the brackets. 
Using the large-A asymptotic expansion of the differential dS, 



where Cl is the contour described above (see also Fig- EJ) ? which by convention encircles the logarithmic 
cut between two infinities on two sheets of the Riemann surface and passes through the last, nth, cut. 

Thus, we proved that dF/dtQ coincides dTo/dto, (|43|) . This completes the proof that the derivatives 
of F and J-q w.r.t. all 5,'s and tk, k = 0,1, ... coincide. Therefore, the planar limit J-q of the 1MM free 
energy is indeed the SWW prepotential or quasiclassical tau- function. Note that this identification of 
the matrix model free energy and the SWW prepotential is crucially based on formula (|72j) which fixes 
solutions to the loop equations (see ss. 1.2-1. 3). Moreover, making higher genera calculations, we shall 
solve the loop equations with similar additional constraints that fix the solution, see s.4 and formula 
(1HT1) below. 

We now introduce the (complete) set of canonical variables {Si, i = l,...,n — 1; to; tk, k = 
1, . . . , ?n}, which we uniformly denote {ti} (in what follows, Latin capitals indicate any quantity from 



dS(X)\ x 



V'(X)dX + -^-dX + 0(A~ 2 )dA 



(76) 



we obtain that (S(b n ))+ just cancels the term V(b n ), and we eventually find that 




(77) 
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Figure 3: Cuts in the A-, or "eigenvalue," plane for the planar limit of 1MM. The eigenvalues are 
supposed to be located "on" the cuts. The distribution of eigenvalues is governed by the period 
integrals Si = § A . p(\)d\ along the corresponding cycles, and the dependence of the free energy on 
"occupation numbers" Si is given by quasiclassical tau-function ^ = § B y(X)d\. We must add the 
logarithmic cut between two copies of the infinity on two sheets of the hyperelliptic Riemann surface 
in order to calculate the derivative w.r.t. the variable to. 

this set). From (|51j) . (|57jl. and (|7.'-{|) . we then obtain the general relation 

ddS _ HjjXjdX 

m7 = dni = ~WT' m 

where Hj(X) are polynomials. 

Asymptotic formulas (|57j) and (|74|) imply that derivatives of all the quantities d£li w.r.t. any 
parameter are regular at infinity and may have singularities only at the ramification points /x a of 
reduced Riemann surface (|21j) . 



2 Second derivatives of free energy 
2.1 Bergmann bidifferential 

In the previous section we demonstrated how to calculate one-point resolvent. This required the 
knowledge of the first derivatives of the matrix model free energy. These could be expressed through 
the local quantities (integrals of differentials) on Riemann surfaces (and allowed us to associate our 
construction with the Whitham system). 

In this section, we apply this procedure to the two-point resolvent, which requires the knowledge of 
the second derivatives. Here, instead of differentials with some prescribed properties of holomorphicity, 
the main object we need is a ^differential: the Bergmann kernel (canonically normalized bidifferential 
in Fay's terminology) which is the bi-differential on a Riemann surface T, g being the double derivative 
of logarithm of the prime form E(P, Q) such that it is symmetrical in its arguments P,Q £ S s and 
has the only singularity at the coinciding arguments where it has the behavior (see |43j ) 

B{p ' Q) = ( (e(p)-e(Q)) 2 + \ Sb{p) + o(1) ) (79) 

in some local coordinate £(P) in the vicinity of a point P £ C; Sb(P) is the Bergmann projective 
connection (Sb(P) transforms as a quadratic differential plus the Schwarzian derivative under an 
arbitrary variable changing; this transformation law is the same as for the energy-momentum tensor 
of the (free) scalar field, see [44] ) . As it stands, we can add to (|79[1 any bilinear combination of 
Abelian 1-differentials duoi\ we fix the normalization claiming vanishing all the integrals over A-cycles 
of B(P,Q): 

I B(P,Q) = 0, for i = l,...,g, (80) 

J A, 
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and, due to the symmetricity property, the integral may be taken over any of the variables P or Q. 

The prime form, a fundamental object on the Riemann surface is defined as follows. Consider 
the Jacobian J, which is a /i-dimensional torus defined by the period map of the curve S ff . Recall 
that the Abel map S g ^ J : P — > x(P) = j Jp dui\, where Pq is a reference point, set into the 
correspondence to each point P of the complex curve the vector in the Jacobian, and we also introduce 
the theta function Oj q j(x) of an odd characteristic [a] that becomes zero at x = 0. Introduce also 
the "normalizing" holomorphic half-differentials h a (P) determined for the points of P £ T, g and 
characteristics a by 

The explicit expression for the prime form E(P, Q) that has a single zero on the Riemann surface T, g 
then reads 

s Qt a ](x(P) -x(Q)) , N 

h a {P)h a {Q) 

while the Bergmann kernel is just 

B(P,Q) = d P d Q log E(P,Q) (82) 

which can be immediately stated by analytical properties and zero ^-periods. 

Note that, up to a holomorphic part, the Bergmann bidifferential is nothing but the scalar Green 
function, see ^Ql H3 HE] • Another useful Green function, that is, the fermionic one, ^ e (P, Q) is defined 
to be a holomorphic 1/2-differential in both variables but the point P = Q where it has the first order 
pole with unit residue. It also depends on the choice of theta-characteristics e (boundary conditions 
for the fermions) and is manifestly given by 

V (P0) = &e(x(P) - X(Q)) 

eK e e (0)E(P,Q) v ; 

The square of the Szego kernels and the bi-differential B{P, Q) are related by the identity [331 (Propo- 
sition 2.12; see also Appendices A,B in [TO} Il2*]): 



* e (P,Q)*_ e (P,Q) = B(P,Q) + du^tkjjiQ)-^- log & e (x) 



(84) 

2 = 



This allows one to express B(P,Q) through the square of the Szego kernel (note that, for the half- 
integer characteristics, —e is equivalent to e). 

As we shall see in the next subsection, the Bergmann kernel generates the differentials (if^. There- 
fore, formula (|82|) would allow one to express these latter through the prime form. Similarly, the bipole 
differential (|57|) can be rewritten through the prime form as 

^-"^W^) (85> 

The primitive of differential (|85|) (which we need in what follows) then obviously develops the loga- 
rithmic cut between the points of two infinities on the Riemann surface. 



2.2 2-point resolvent 

Now one can easily express the 2-point resolvent Wq(\, (i) in terms of B(P,Q) on hyperelliptic curve 
(|21|) . where we now use the hyperelliptic co-ordinate A = £(-P) and fj, = £(Q). Indeed, let us use ©, 
(|62|) and (JHH) to obtain 

W (A, fi)dfid\ = k^xm Tes oo+x k dni(x) = -r^dflifji) (86) 

k,i>o V l>0 
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where d£lk = dflk — \dx k is the meromorphic second-kind Abelian differential with the only singularity 
at oo_, where it behaves like 

dCl k (fi) = -k U k ~ l + O(l)) dfi, for /i -► oo_, k > (87) 

and has vanishing A-periods £4. dVL^ = 0. Therefore, Wq(X, fi) is holomorphic everywhere if both A 
and /x correspond to the points on the same sheet, but it develops the second order pole at // = A, 
where two points are located on different sheets. Indeed, taking into account the only non-holomorphic 
part gives 

k\ k ~ X 1 

Wb(A )A i)~ -J2 



Besides, the evident normalizing condition, fixing the holomorphic part, immediately follows from 

I W {\, fi)dfi = <f W (X,fi)dX = (89) 

J A, J Ai 

due to (|T2|) and since Wq(X, n) is symmetric in A and n by definition. 
Therefore, we finally come to the formula for the 2-point resolvent, 

W (X,fi)dXd» = ^^-dXdfi = -B(P,Q*), (90) 

where we have introduced the *-involution between the two sheets of the hyperelliptic curve C, so that 
Q* denotes the image of Q under this involution. The only singularity of Q9U|I . for a fixed point P on 
a physical sheet, is at the point Q — > P* on the unphysical sheet with fJ,(Q) = X(P*) = X(P), while 
on the other sheet it is cancelled under change of sign of y. 

Therefore, in order to calculate the 2-point resolvent, one needs to write down the Bergmann 
bidifferential on the hyperelliptic curve manifestly. In principle, it has several different representations 
(one of the most hard for any further treatment is given by formula (5.20) in ^U], borrowed from |46|). 
The simplest one can be obtained using formula (|84[). Indeed, a simple hyperelliptic representation of 
the Szego kernel exists for the even non-singular half-integer characteristics. Such characteristics are 
in one-to-one correspondence with the partitions of the set of all the 2g + 2 ramification points into 
two equal subsets, and {//"}, a = 1, . . . , g + 1, y±{X) = U 9 a t\(.^ ~ Mo)> le - v( X ) = V+Wv-W- 

Given these two sets, one can define U e (X) = y y + £j . In terms of these functions, the Szego kernel is 
equal to jSl EZj 

2y/U e {X)U e (n) X-fi 

Square of this expression, due to (|84|). leads to the manifest expression for the singular part of 2-point 
resolvent or the Bergmann kernel, i.e. 

Wo(X, u)dXdu = )y~(^) + y+(/J-)y-( ) V V ^ ^_ _|_ holomorphic part (92) 

4y(A)y(/i) (A -m) 2 

where we choose the sign in front of 2y{X)y{y) in the numerator so that the pole lies on the unphysical 
sheet. The holomorphic part is fixed now by the condition of zero ^4-periods. 

However, in our further calculations we need another, completely different expression for the 
Bergmann kernel [HI]. It can be most immediately obtained from the loop equation. The loop 
equation for the 2-point resolvent has the form (see, e.g., formula (1.3.40) in |20| ) 

V'(X)Wo(X,») - r v (X)W (^ = 2W (\)W Q (\,n) + ^ W °W~ W <>M (93) 
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where the operator ry(A) is defined in Q17|). i.e. on a hyperelliptic curve y 2 = R(x), e.g. (fT^|) . one gets 



W (A, M ) 



y(A) 



5 Wq (A) - W (ji) 



y(A) 

2(A-^ 2 ' 2y( M ) 



+ 



1 



A — /i 

1 dlogi?(A) 1 



+ MA)W (m) 

»V(M)logi2(A) 



(94) 



SA 



.(A -a*) 2 2(A-/i) 

One can check by straightforward calculation that this formula leads to a symmetric expression (see 
(III.2.6) in [20]) 



Wq(X, fi)dXdfj, 



V'(jj,)V'(\) + \ (AQ(m) + M<3(A)) + c - y(^)y(A) dXd^i 



+ holomorphic part (95) 



2y{X)y{ l i) (A - ,u) 2 

which is a particular case of (|92[). when parameterizing the hyperelliptic curve as 

y 2 (A) = y /2 (A) + AQ(A) + c (96) 
In the case of degenerate Riemann surface Q21JI . one obtains (see [HI]) similarly to (|94|) 



Wq (A, n)d\dfi 



dXdfi 



+ 



2n r 



a=l 



(/X - A)(A - H a ) \~Ha 



+ 



dAd//, (97) 



where C a ((A) = Y^=® £a,iH l are polynomials in /i. They can be unambiguously fixed by the require- 
ments of absence of the first-order poles at A = fi and zero ^4-periods. 



2.3 Calculating C a (fi) 

Although one has now the explicit formula for the 2-point resolvent (=Bergmann bidifferential) on 
the hyperelliptic surface in terms of branching points (|97|). it contains the polynomials J~. a (/j,) defined 
by the implicit requirements. 

To find effective formulas for these polynomials, let us set A = fi a in (|97|). and introduce the 
notation 

y*W= IHix-np), y a = JH^a-vp) (98) 

V V 

We will also denote by square brackets the fixed argument of the Bergmann bi-differential (in some 
local coordinate), and consider B(P, \p]) as a 1-differential on C. From (|9U() . one has 

- 2B(P, K]) = \^vld\ + \ V hrm dX (") 
2 y(X) 2 (A - fi a )y{\) 

From (JHUJ), we have 

/ SMit = - 1 » (ioo) 

and integrand in the l.h.s. is a linear combination of canonical holomorphic differentials duoi (|51|). so 
that 

' dA 



^ .M< (A-// Q )y(A) 

and, in particular, 

n— 1 „ j\ 

^ /i, (A - Ma)y(A) 
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Another equivalent representation we will need in sect. 0J is 

n— 2 n— 1 



£4*,tf4 = -E# n \~(U dX ' <* = l,...,2n. (103) 



Indeed, let us introduce the quantities 



f A* -1 
J a, y{X) 



Then, for the canonical polynomials -fffc(A) = J^^i Hi^X 1 1 , k = 1, . . . ,n — 1, related to the canoni- 

n-l 



cally normalized differentials (|51|) . i.e. <iA = one obviously has 



E ^'.J-Hj,* = fy,* for j, A; = 1, . . . , n - 1. (105) 



Therefore, for all > such that j — 1 — > 0, 



Then, 



n-l „ \j-k—l 

E^,.£ w ,A = (106) 



^ 4 (A-^)y(A) dA " 

y(A) 

and, because of ()102|) . we finally arrive at ()103|) . The above formulas mean that for the Bergmann 
kernel on hyperelliptic curve y 2 = R(x) one can write 



2/i 



dXd^ y(Xl / 1 1 v 



2(A-//) 2 2y(/i) — A) 2 2 



Q = l 



Wo(A, n)dXd/i = 

(/,-A)(A-^) " ^ ^ £ (x- M «) 2 l/(x) 



dXdfi 



dXdu ( R(X) \ dXdfi R'(X) 1 ... ^^-4 , , . / , inw , 
" I " ) ~2(A - fi) y(X)y(p)~2 £ g **<"> A, (, - (l ° 7> 

A particular case of this formula at A = was used in [IS] for solving the quasiclassical Bethe anzatz 
equations in the context of AdS/CFT correspondence. The last term in the r.h.s. of (|!07j) is explicit 
form for the action of the operator ((17)) in the case of smooth Riemann surface, restricted by the 
vanishing period's constraints. 

2.4 Mixed second derivatives 

Another set of relations follows from the general properties of the Bergmann kernel, and can be also 
derived directly from the formulas of sect. ^ To this end, we apply the mixed derivatives d/dV(fj,) 
and d/dSi to the planar limit of the free energy Tq. On one hand, d^o/dSi = § B , dS and using that 

dS{X) = y{X)dX = (V'(X) - 2W Q (X))dX, = -jj^ and formula (gQ), one obtains that 



= / 2B(P, [/*]). 
JBi 
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On the other hand, acting by derivatives in the opposite order, one first obtains dJ : o/dV(fj,) = Wq(/u) = 
V'(fi) —y(fJ') and then d(V'(fi) —y(fj,))/dSi = 2diOi([/j]), or, in the coordinate-free notation, one of the 
Fay identities 

/ B(P,Q) = diOi(Q) (108) 



This means that 



ddS(n) 



WS(A) ■ 



(109) 



where the 6o£/i integrals are taken over the variable A. Now, as dS(fi) = y([i)dX is the generating 
function for the variables £ a = {M, {Aj}, {^ a }, M^}, given by (HU), (El, and ^ (M$ are just ith 
order derivatives of M m _ n (fi) at fi = fj, a ) giving their dependence on Si and tk, one concludes that 
similar relation for the first derivatives holds for each of these variables. Indeed, multiplying (J109|) by 
one then can bring \i successively to Ha's, Aj's and oo to pick up pole terms with different £ a and 
prove that 

9ia I d^a ^ (uo) 



dSi J Bl 0V(X) 

As a consequence, any function G of £ a would naively satisfy the same relation. However, there is a 
subtle point here: in formula (|108|) one could integrate the both sides over an A,-cycle to obtain at 
the r.h.s. 6ij, while the A-periods of the Bergmann kernel are zero. It means that one should carefully 
permute the integrations, since Ai— and .Bj-cycles intersect. Therefore, one should carefully take into 
account the contributions of the intersection points (Ai- and -Bj-cycles for different i and j do not 
intersect and integrations can be exchanged, which perfectly match the Kronecker symbol obtained) 
which does not vanish due to the double pole of the Bergmann kernel. 

In particular, we formally have 

dM(\) _ r (dM{\) d 1 \ 

where the second term in the brackets vanishes unless we have an outer integration over the cycle Ai 
w.r.t. the variable A in the both sides. 

Similarly, one should take care when omitting the second term in the second equality in (|108j) . 
Indeed, suppose we consider Si as a function of £ a . Then, 

dSi(Sa) , i OSi 



*Lwk (112> 



dSj ' J Bi dV(X) 

Here again one should note that Si is the integral of dS over the Ai cycle and take care when exchanging 
integrations. While doing this, the contributions from the omitted term in (|108j) are to be taken into 
account. Note, however, that one can never express Si as a function of a finite number of "local" 
variables £ a . 

Therefore, if one considers the functions G that depends on only finite number of "local" variables 
£ a (the branching points, Mq' — the moments of the model, and, possibly, zeros of the polynomial 
M(x)), the relation 

dG f dG „ . , 

-dX (113) 



dSi J Bl dV(X) 

holds, while, for A-cycle integrals over the Riemann surface (note that the -B-cycle integrals do not 
meet such a problem), one should add more terms in the r.h.s. in order to take into account the 
second, pole term in (|1U8[) . <|111|> . 

Relation (|113[) is therefore valid in all orders of 1 /^-expansion for the 1MM free energy because 
any higher-genus contribution is a function only of fi a and of a finite number of higher moments . 
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Equation (|U7|) implies, for A ~ /x Q ~ fi, in local coordinates d\/\/X — \x a = 2d(y/X — and 
dfj,/^fi — fx a = 2d (^/ n — fj, a ) the relation 

£a(Ha) = B([n a ], [n a ])\nonsing. = S B (Ha)- (114) 

Applying now (|108j) and (|103|) . one immediately comes to 

CM = SbM = £ / / (us) 

r^JAiJBi A-/i a 

in the local coordinates associated with the hyperelliptic Riemann surface l|21j). 

3 WDVV equations 

The general form of the Witten-DijkgraafoVerlinde-Verlinde (WDVV) equations |24ll25j is the systems 
of algebraic equations 

TiT- ] 1 Tk = TkT- ] 1 T u V I,J,K (116) 

on the third derivatives 

- d£k s ^ k (nT> 

of some function These equations often admit an interpretation as associativity relations in 

some algebra (of polynomials, differentials etc.) and are relevant for describing topological theories. 

As WDVV systems are often closely related to Whitham systems, a natural question is whether 
the corresponding 1MM free-energy function satisfy the WDVV equations? This was proved in |31j . 
where it was shown that the multicut solution the 1MM satisfies the WDVV equations as a function of 
canonical variables identified with the periods and residues of the generating meromorphic one-form 
dS P|. The method to prove it consists of two steps. The first, most difficult, step is to find the 
residue formula for the third derivatives (|117[) of the 1MM free energy. Then, using an associativity, 
one immediately proves that the free energy of multi-support solution satisfies the WDVV equations if 
the number of independent variables is equal to the number of branching (critical) points in the residue 
formula. We show here that the statement holds in the case of arbitrary potentials for a fixed-genus 
reduced Riemann surface. 

In sect. 13. ll we derive the residue formula for the third derivatives of the quasiclassical tau- function 
for the variables (the generalized times) tj associated with both the periods Si and residues t/. of the 
generating differential dS. In sect. 13. 2| we prove that the free energy of the multi-interval-support 
solution Toiti) solves WDVV equations (|116|) as a function of the subset {t a } C {ti}; the total number 
of t a must be fixed to be equal to the number of branching points in the residue formula for the third 
derivatives (|117|) in order to make the set of the WDVV equations nontrivial. 



3.1 Residue formula 

Because all the quantities dfti ((75)) depend entirely on the reduced hyperelliptic Riemann surface (|2"Tj) , 
their derivatives w.r.t. any parameter must be expressed through the derivatives w.r.t. the positions 
of the branching points fj, a . So, calculating derivatives w.r.t. fi a is the basic ingredient. Note that 
although the differentials are regular at the points fi a in the local coordinate dy ~ d\/\/\ — /j, a , the 
derivatives ddQi / d\i a obviously develop singularities at A = fj, a , and we must bypass these singularities 
when choosing the integration contour as in Fig. |5J 

Let us now derive the formulas for the third derivatives d 3 J r o/(dtjdt jd^ a ) = Tija-, following [§ll3~T]. 
Consider, first, the case where the "times" ij and tj are Si and Sj, and = Tu = Tij. We note 

that the derivatives of the elements of period matrix can be expressed through the integral over the 
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"boundary" dTi g of the cut Riemann surface T, g (see Fig. |2J), Indeed, because of the normalization 
condition § Ai dujj = Sij, we have § A d a du>j = 0, so that 



(p d a du>i = I (p duij <z> d a du>i — <p doj~ <p d a du>i 
Jb j l^X ^ Ja i Jb i Jb i Ja i 

^2 ( / u+d a duJi - / ujjd a duji J - 

l— l \J Bi JBi / 



1=1 

ujjd a du>i (118) 

where ojj = J dujj are Abelian integrals and we let lo^ denote their values on two copies of cycles on 
the cut Riemann surface in Fig. [21 

We must now choose the cycles A[ and B\ bypassing all possible singularities of the integrand (in 
this case, the ramification points ug). Expression (|118f) can be then evaluated through the residue 
formula 

d a Tij = - ujjd a duji = ^res d A=o (d a LOjdLOi) . (119) 



The proof of this formula for generic variation of moduli can be found in |31j . Here we will adjust it 
for the class of variations in terms of the branch points {fi a }, i.e. to the class of hyperelliptic Riemann 
surfaces. 

Before evaluating this sum of residues, let us consider the case of meromorphic differentials. Then, 
using formulas (|62|) and (|77j) . one obtains 

' ((n fc )+.odnO , M>o, (120) 



dt k dti 2 Jc L 

where (flk)+,o is the singular part of fifc at infinity, i.e., it is X k for k > and the logarithmic function 
for k = 0. Because d(fLk)+,o/dfJ> a = for k > 0, we have ddVt k /d^ a = d(dflk)- /dfi a , where (dQ^)- 
is the holomorphic part of dVtk at infinity (and the expression (f^)- is therefore meromorphic for all 
k > 0). We then have 

d 1 



.. )? ((n k ) +:0 dn l ) = l<f ((n Jfc ) +i0 a a (dn { )_) = 

dua 2 J Cl 2 J Cl 

= -- I {dtt k d a (Qi)) = -resoo {d£l k d a tti) . (121) 
2 Jc L 



The last expression can be rewritten as 



d ^ lk 'djj~) = Jot, \ tk d~) + ^ r6SM/3 ( d ^ lkda ^ l " > = 

Q (3 = 1 

2n 

= res^ (dn^i) (122) 

(3=1 

because §q S (dflkd a ^li) = following the same arguments as in formula (|118|) and due to normalization 
conditions (1721). 



Further computations for both holomorphic and meromorphic differentials coincide as we need 
only their local behavior at the vicinity of a point A = /ig. Using explicit expression (|78|) . we have 



Hj(/J. a 



8 a nj= (A-/i a )- 1 / 2 + Q( v ^7^), (123) 

li-r^a y^a — 
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for (5 = a, and d a £lj ~ ^A — Hp otherwise. Together with (JTHJ), this means that the only point to 
evaluate the residue is A = Ha at which we have (cf. [33]) 

Sojj a = res^di lid a ilj) = ^ -. r. (124) 

Llp^aKHa - HP) 

Completing the calculation of the third derivative needs just inverting the dependence on the 
ramification points therefore finding d\x a jdtK- Differentiating expressions ()19|) . (|21|) w.r.t. tx for 
computing ((ST?]) we obtain 

ddS _ H K (X)dX _ 

l M /nT ^( A ) i d M m-n(A) 

= -M m _ n (A) }^ — -— dX + — y{X)dX. 125) 

The derivative of the polynomial M m _ n (A) is obviously polynomial and regular at A = Ha- Multiplying 
()125|) by y/X — fi a and setting A = Ha, we immediately obtain 

dfM = H K (fj. a ) ^ 12g ^ 

dt K M m _ n (Ha)Ylp^ a {Ha- Hp)' 

Combining this with (|124j) . we come to the desired residue formula for the third derivative w.r.t. the 
canonical variables tf. 

d 3 ^ ^ Hi([i a )Hj(ijL a )HK(pa) sr^ d^jdQ,jdQ,K dQidQjdQx 
dttftjdtx ^ M m _ n (fi a )l\p^ a (fi a - (ip) 2 ^ reS>la dXdy dXdy 

(127) 



3.2 Proof of WDVV equations 

Given residue formula (|127|) . the proof of WDVV equations (|116j) can be done, following |26|-|29|. 
by checking associativity of the algebra of differentials d£li with multiplication modulo y(X)dX. This 
algebra is reduced to the algebra of polynomials Hj(X) with multiplication modulo y 2 (A) which is 
correctly defined and associative. The basis of the algebra of Hj(X) obviously has dimension 2n and 
is given, e.g., by monomials of the corresponding degrees 0, 1, . . . , 2n — 2, 2n — 1. The study of such 
algebras can be performed even for non-hyperelliptic curves, the details can be found in |28| I49j . 

Another proof is even more simple and reduces to solving the system of linear equations |5L)l I3L)| . 
To this end, we first define 

Hi (Ha) 



so that (|127|) can be rewritten as 



(128) 

Mrrl-niHa) Up^Ha ~ Hp) 2/3 

T^ I jK = y j 4>Wj4>' X K (129) 



Now let us fix some index Y and consider the following multiplication 



Y. C u )K <t>K<t>Y, Va (130) 

K 



the structure constants C^p K being independent of a. One can equally look at this as at a system 
of linear equations for Cfj at fixed values of / and J. If this system has a solution, ()13U|) gives rise to 
an associative ring, with the structure constants Cfj satisfying (associativity condition) 

Kt( c 'P)>( c fT L ( c nlr (131) 
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Now, the solution to (|13Uj) is 

c£P* = £*?*j(*fo*) (132) 

a 

where is the matrix 18 inverse to (p^cpy- This solution exists if the number of vectors (variables 
K) is greater or equal the number 2n of the branching points u a . 

The other important condition is the invertibility of the matrix </>^-</>y, or the matrix cjfy (we 
suppose 4> Y 7^ 0) which ensures the nontriviality of the WDVV relations. For this, we must require 
the number of vectors (pi to be less or equal the number In of their components. We therefore obtain 
the following two conditions |3()| : 



the "matching" condition 

#(/) = #(a); (133) 

and 

the nondegeneracy of the matrix <pj (see the proof in sect. I4.3JI : 

detail #0 (134) 

la 



Now, using (|13()|) . one rewrites (|129j) 

Fqjjk = £ Cfj^^m = E (135) 

coming to the matrix formula that express the structure constants through the third derivatives of 
the planar limit free energy 

C ( P = TiT Y x (136) 

where we denoted Tj the matrix with element JK equal to !Fq,ijk- Now substituting (|136|) into (|131j) . 
one immediately arrives at ()116|) . 

Thus, we established that conditions (jllfi[) require the number of varying parameters {£/} to 
satisfy matching condition ()133|) . We let {t a } denote these "primary" variables. It is convenient to 
set classical "primary" variables w.r.t. which WDVV equations (|llfi|) hold true, to be the parameters 
Si, i = 1, . . . , n — 1 = g, to, and t/. with k = 1, . . . , n keeping all other times frozen. 

Below, in sect. 14. 3| we interpret the answer in genus one in terms of the determinant relation for 
the third derivatives Tuk- 



4 Higher genus contributions 

The solution M^i(A) to the loop equations in the multicut case was first found by Akemann [HI]. 19 He 
also managed to integrate them in to obtain the free energy T\ in the two-cut case. The genus-one 
partition function in the generic multi-cut case was proposed in [^E2], where it was observed that 
the Akemann formula coincides with the correlator of twist fields, computed by Al.Zamolodchikov 
|53| . This produces cuts on complex plane and gives rise to a hyperelliptic Riemann surface, following 
the ideology of jl^j, some corrections to this construction are due to the star operators, introduced 
in [21]. In this section, we present (see also |55j ) the derivation of the genus-one correction based on 
solving the loop equation, and generalizing Akemann's result for the partition function to arbitrary 
number of cuts. 

18 We consider it as a matrix of indices K and a, while the reference index Y is implied as a silent parameter of the 
whole consideration. 

19 The universal critical behavior of the corresponding correlation functions was discussed in 0. 
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4.1 The iterative procedure 



Iterative solving the loop equations. Thus, now we are going to determine higher genus contri- 
butions. We do this iteratively by inverting the genus expanded loop equation Our strategy will 
be to construct an integral operator dQ inverse to the integral operator K — 2Wo(A). 

Acting with this operator onto the both sides of the loop equation ea. ((TT|) . one recursively produces 
W^(A) for all genera like all the multi-point resolvents of the same genus can be obtained from W^(A) 
merely applying the loop insertion operator gy^ ■ 

However, there is a subtlety: the operator K — 2Wq(X) has zero modes and is not invertible. There- 
fore, solution to the loop equation is determined up to an arbitrary combination of these zero modes. 
Hence, the kernel of the operator K — 2 Wo (A) is spanned exactly by holomorphic one-differentials on 
the Riemann surface (|21|). 

In order to fix this freedom, we assume W^(A) is expressed exclusively in terms of derivatives Jy^ 

and , which, as we show in the next paragraph, fixes a solution to the loop equation. It is a 

natural extension of the normalizing property (|6U|) to higher genera and can be ultimately written in 
the form 

' 11 dX = <b W h (X)d\ = Vi and for h>l. (137) 



T Ai dV(X) Ja 
Now we claim that the integral operator 

»(/)(*)=/ ^^-^-/M (138) 
Jc„ a 2m dX y(p) 

is an inverse for the operator K — 2 Wo (A) in the space of rational functions /(/i) with poles at the 
points fj, a only. 

Here the one-differential dQ(X,/j,) w.r.t. the first argument A 20 is the primitive of the Bergmann 
kernel B(X,n) w.r.t. the argument /i. Obviously, it is a single-valued differential of A with zero A- 
periods on the reduced Riemann surface and is multiple-valued function of which undergoes jumps 
equal to dwi(X) when the variable \i passes through the cycle Bi (cf. with (|l()7jl ): 

{X-n)y(X) y(A) J Al (C - Nv{0 

Contours of integration over the cycles Ai must lie outside the contour of integration C^ a encircling 
the branching point \i a in (|150|) . 

Moreover, this operator obeys the property 

^(/)(A)dA = (140) 

and, therefore, respects condition (|137[) . Therefore, one has to solve the loop equations inverting 
K - 2W (A) exactly with dQ. 

Then, the calculation immediately validates the diagrammatic technique |32| for evaluating mul- 
tipoint resolvents in 1MM. Indeed, representing dQ(X,fJ,) (for A > /i) as the arrowed propagator, the 
three-point vertex as dot in which we assume the integration over fj,: • = § ^^y, we can graphically 
write solution to (jllj) since 



W h (X) = dQ 



h-l 

E WA'(-)Wft- fc /(o + Wfc-i( 

Lh'=l 



(A). (141) 



It is the function dS(X, fi) in the notation of 
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Then, representing multiresolvent Wh'(Xi, . . . , A&) as the block with k external legs and with the index 
h', one obtains 




(142) 

which is just the basis relation for the diagrammatic representation (|141|) . Here, by convention, all 
integration contours for the variables /i lie inside each other in the order, established by arrowed 
propagators dQ(fii, The other, nonarrowed propagators are W(X,fi) = ^-dQ(X, fi). All ^4-cycle 
contours of integration in dQ(X,n) are outside the contours of internal integrations over /ij- variables. 



Choosing a specific basis. In order to prove the claims of the previous paragraph, first of all, we 
change variables from coupling constants to special moment functions which allows one to apply higher 
genus machinery nonperturbatively in coupling constants tj. This machinery turns out to involve only 
on a finite number of the moments (JUJ). (Recall that (f23 > (f!H |l implies M^P = M{ji a ).) 

Now let us fix the generic analytic structures of the 1- and 2-point resolvents. First of all, note 
that Wo(X, in (|94|) is invariant w.r.t. the involution y — * —y that permutes physical and unphysical 
sheets. Moreover, Wo(A,A) is a fractional rational function of A with poles at the points fi a only. 
Further, look at the loop equation, (|11|) and rewrite it as 

h-l o 

y{\)W h {\) = [V(X)W h (X)] + + W h ,(\)W h - h ,(\) + ^yW h _!(A), (143) 

It follows from this formula that the 1-point resolvent Wi(A) is also fractional rational function of A 
with poles at the points fj, a only divided by y(X), i.e. it changes sign under permuting the physical 
and unphysical sheets. This procedure can be iterated with the loop equation written in the form j5U] 

h 

V'(\)W h (\,\ 1 ,...,\ n )=r v WW h (\i,...,\n)+Yl E W h '(X,X 1 ,...,X ni )W h - h >(X,X 1 ,...,X n2 )+ 

h'=0 ni+n2=n— 1 

i 9 ^(A,Ai,...,Aj,...,A n ) - VF fe (Ai,...,A n ) d 

+ \ d\ t — t + wm Wh - i{x > Xi > Xn) (144) 

In particular, using the known analytic structure of Wo(A,//) one easily checks that the 3-point resol- 
vent Wq(X, fj,, v) is odd w.r.t. to the involution y — * —y and then, with the knowledge of structure of 
Wi(X) and Wo(A,//), one can use the loop equation 1)144(1 to prove that W\(X, fi) is even w.r.t. the 
involution etc. The final result is that all (2n)-point resolvents Wh are even, while all (2n + l)-point 
resolvents Wh but Wo(X) are odd w.r.t. the involution. This means that all Wh(X) but Wo(A) are 
fractional rational functions of Aj with poles at the points fj, a only divided by y(X). Moreover, the 
r.h.s. of eq. (|11|) is similarly a fractional rational function of A having poles at \x a only, i.e. one 
should naturally choose a specific basis \a (A) defined by the property that, for the integral operator 
in ea.lfTTj). 

(K-2W (X))x { a h \X)= 1 

[A — fj, a ) 

lb = 1,2,..., a = l,...,2n. (145) 
Then, ^^(A) must have the structure 



3/i-l 2n 

W h (X) = E E <ixL fc) (A), h > 1, (146) 

k=l a=l 
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(k) (k) 

where A a h are certain functions of /xg and the moments . As the order of the highest singularity 
term 1/((A — // Q ) 3 ' l-1 y(A)) in Wh(X) is insensitive to a multi-cut structure 21 , Wh(X) will depend on at 
most 2n(3/i — 2) moments, just like the one-cut solution case |35j . 

One could define a set of basis functions Xa\X) recurrently, as in [31], however, here we 
present another technique inspired by [SI]. We first calculate the quantities d yfy and dv fy ■ 

Using the identity qv^(a) ^(m) = ~ (A-//) 2 anc ^ re P resen tation (|24[). one easily obtains 
^ - (jfc + 1/2) f 1} ^ 1 



aV(A) 7 y V Q 5F(A) (A - M«) fc+1 i/(A) 

+ - y y - ( - 

+ 2^( / , /3 - / x Q )H ^ SU(A) (A-^)y(A)J 1 ' j 

a = 1, . . . , 2n , A; = 1, 2, ... . 
Note that the general structure of this formula is 

om^ = dfi a v d / i \ 

ay (A) W dF(A)^ W dF(A) <9A V (A - /%) fc y(A) J ' 1 j 

In order to calculate the derivative of the branching point /x a w.r.t. the potential, one can note 
that Wq(X, fi) = \ d ^ V qv^x)^ and bring the variable [i in this expression to fi a . Then, from (|U7|) and 
(|1U1|) . one obtains that 

a 9V(A) (A-/v>y(A) ^ y(A) 7* (£ - M a)y(£) 1 j 

It immediately follows from these formulas and (|1U0|> . (|101|) that integrals over ^-cycles of both 

o AA/fC'') _____ 

and gyfh vanish and one, therefore, arrives at (|137[) . 

Using the above conditions and formula Q137|) . we can now invert the operator K — 2 Wo (A) when 
acting on basis monomials (A — fJ. a )~ k . That is, we are going to check that the basis X^(A) vectors 
(|145|) are generated from these basis monomials by the operator dQ 

JC„ a 2vrz y(fi) dX (p - fi a ) k - V > 

First few basis functions are easy to obtain from (|149j) and (|139[1 : 

vWm = d ^ a a = l 2n 

«(2)m = -- - loglM^I- 
Xa W 3 dV(X) 1 

1 2n 5 

-3E^y lo gK-M ( 151 ) 



'This can be also stated from the analysis of the loop equations as above. 
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Proof of (|138|) . First, let us demonstrate that the action of the operator dQ defined in (|15U|) on 
any basis function (A — fJ> a )~ k inverts the action of K — 2Wq{X) up to the zero mode content. For this, 
let us consider the expression 



(K - 2W (X)) 



d/i dQ{\n) 



a I 



dw V'(w) f dfj, dQ(w,/j,) 1 



c Vw 2m A - w Jcv^ 2m dw - fi a ) k 

-2W (X)i 1 



p, 



2m dX y{n){ji- /i a ) k ' 



Taking into account that the contour ordering is such that Cd w > Cx> M m the sense that one lies inside 
the other and evaluating the integral over w by taking residues at the points w = X and w = oo, we find 
that the result of the residue at w = X combines with the second term to produce V'(A)— Wo(A) = y(X) 
while we can replace V'(w) by y(w) when evaluating the residue at infinity due to the asymptotic 
conditions. That is, we obtain 

/ du, dQ(X,u) 1 / dw y(w) f dfi dQ(w,a) 1 

9. oZ7 7\ 77377 TTTfc + 



Ci , m 2m dX y(fi)([i - [i a ) k Jc^ 2m A - w Jc^ 2m dw y(n){n - ^ a ) k 
dw y(w) f dfj, dQ(w,fj,) 1 



t, w 2mX-w J Ct , u 2m dw y{n)(n - ^ a ) k '' 



where the contour ordering is such that A > Cx> w > Cx> M • We now want to push the integration contour 
for [i through the integration contour for w. After it, the obtained integral over [i vanishes as the 
integrand is then analytic everywhere outside C?> . Thus, contributions come only from the pole at 
w = fj, of dQ(w,fi), which contributes when pushing the contour Cx> M through Cx> w , and from the 
multiple-valuedness of dQ{w,fi) w.r.t. the variable fi. Note, however, that these latter contributions 
are always proportional to duji(w) = — ^j^y— ■ That is, we have 

dw y(w) 1 
2m X — w y(w)(w — /J, a ) k 

dw d/i y(w) Hi(w) 1 



+ const 



A>Cx, u >Cp„ 2ni 2m X - w y(w) y{p)(n - n a ) k ' 



where the point A lies outside the integration contour in the first term and the integral over w in the 
last term vanishes because the integrand 

y{w)Hi{w) M{w)H l {w) 



(A — w)y(w) X — w 

is obviously regular everywhere inside the contour Cx> w (recall that both A and fi are now outside this 
contour). Upon integration, the first term obviously produces (A — /-i a )~ k , which completes the first 
part of the proof. 

Next, note that conditions (|137|) hold automatically for any function /(A), having singularities 
only at /i a , transformed by the operator dG, (|14()|) due to the normalization properties of the kernel 
dQ(X,fi) H139j) . Therefore, the result of the action of this operator can be always presented as the 

o o * *- ( fa ) ^^^^^^^^ 

linear combination of the functions gyjjj and q V "^ , which completes the proof of formula (|15fl|) . 
4.2 Calculations in genus one 

Now we invert the loop equations for genus h = 1 and integrate them to obtain the genus one free 
energy. We need, in this case, only Xa\x) and x« (A) (see (|146j0 which we already have, (|151[) . and 
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ea. lfTT|) reads 



(K - 2W (X))W 1 (X) 



a 



Given Wo(A) (|26]> . the r.h.s. becomes 
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(152) 
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Here we took into account that regular parts coming from ^ l , Z = 1, . . . ,n— 2, vanish for Wo(A, A) = 

8V{\) Wo (A) to satisfy the correct asymptotic behavior, and we can just replace A' by fi l a in numerators 
of such expressions. The result for the one-point resolvent of genus one with n cuts can now be easily 
obtained using Eq. (|151|) . 
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Now one should integrate Q154j) in order to obtain ^i. While integrating the first two terms in the 
r.h.s. is straightforward, the term with the zero modes requires some more work. Using formulas (jlU3j) 
and Q1U5|) . we see that the last term in (|154f) is just 



\2K d 



log det (To j 

i,j=l,...,n-l 



av(xy 



i.e., 



T X = -^log (pW • A 4 • (.. -le-,, ; -,/)' 2 ) 



(155) 



where A = Yi\<a<p<2n{Ha ~ Hp) 1S the Vandermonde determinant. This is our final answer for the 
genus-one partition function 22 . 

22 One has to compare this answer with the formula proposed in |56l 1241 157| for the one- loop (toric) corrections in 
topological theories, 



Fi(ti) = — log det 



dtadtpdt 



-d x t~, 



+ G({t a }). 
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4.3 Genus one free energy and determinant representation 



Let us now discuss the expression for the genus one free energy. First of all, notice that (|155|) reproduces 
the calculation of 34j for the two-cut solution, up to a modular transformation permuting A— with B— 
cycles. This should be a surprise, since we put throughout our calculation in sect. 14.21 the constraint 
that A-periods ()3U|) of the generating differential (|67[) are constant under the action of the operator 
, see (|72|) ■ On the contrary, in [33] Akemann imposed the condition of vanishing B-periods of dS, 
corresponding to equal "levels" IT in different wells of the potential j^H] or additional minimization 
of the free energy (|39[) at the saddle point w.r.t. the occupation numbers ([41 j) . In fact, since neither 
the answer (j!55j) nor intermediate calculations contain any manifest dependence on particular values 
of the periods, one can equally put all B-periods (during the calculation of [33]) fixed to be arbitrary 
non-vanishing constants and, therefore, come to a modular transformed counterpart of our choice of 
the normalizing cycles. 

Under condition of constant B-periods (|44 |) -(|47 |) . as it was stressed in [23], the matrix Gi j of the 

j4-periods of is replaced by the matrix of the corresponding B-periods. This is the only difference 
y\ x ) 

with the result of [S3]; certainly formula <|155|) reproduces the answer of |51M52j 23 for generic multi-cut 
solution. 

The fact, that the only result of interchanging A- and B-cycles is the interchanging of the corre- 
sponding periods in det a implies that e^ 1 is a density and not a scalar function on moduli space of 
the curves. Indeed, when exchanging A- and B-cycles, detcrjj is multiplied by detrjj - the determi- 
nant of the period matrix of the curve. In order to compensate this factor, e^ 1 must be transformed 
under such transformation with the additional factor (detrjj) 1 ^ 2 , as follows from (|155|) . Then one 
immediately comes to the above observation: exchange of A- and B-cycles results only in replacing 
the corresponding matrix Gi j in (|155|) . Note that such behavior of e^ 1 indicates that it is a section 
of determinant bundle DET9 over the moduli space, where the <9-operator acts on the sections of a 
non-trivial bundle on a complex curve of matrix model. One can find that determinant det' dj of the 
<9-operator (with some fixed basis of the zero modes), acting on j-differentials, is proportional to det a 
for j = 0, 1 but for other values of j it typically does not contain the factor deter, still transforming 
non-trivially under exchange of A- and B-cycles. It was proposed in |52| that, in order to match 
the proper behavior under modular transformations, the operator dj should act on twisted bosons on 
hyperelliptic curves, then e^ 1 actually equals to its determinant. Besides, one also needs to add some 
corrections from the star operators [52| I54| that do not contain det a factors and cannot be restored 
by modular covariance of the answer; these are necessary to obtain the correct result (|155[) . 

There is another important point that differs between our formula and the result of [33] • Namely, 
while in [34| it was possible to add any constant to the final result, not spoiling the solution to the 
loop equation, we can add to (|155j) an arbitrary function of occupation numbers Si, since in contrast 
to |34j with no free parameters, we keep Si's arbitrary. 

One can partially fix this arbitrary function in the free energy by imposing requirement of smooth 
behavior of T\ under degenerations of the surface. To this end, let us shrink one of the cuts, e.g. bring 
^2 to Setting fj,2 — m = e — > 0, we can easily check that 



Jn) 1 , 



In m~n 
a=3 i=l 



+ f{ n - 1] + 0(e) (156) 



In order to compensate the first term that spoils the smooth degeneration, one suffices to add 

n n—1 

II log S n = t -Y,Si (157) 



12 . . 

t=i t=i 



Note that this term is out of control in the conformal field theory approach of [511 152j . On the other 
hand, it could be also compared with the matrix model calculations of [S1H3]- An arbitrary function 

23 If restoring in |51| the determinant term det a, omitted from the answer. 
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of Si comes there from different normalizations of the matrix integral. In particular, the normalization 
in [SJ[23] corresponds just to (|155|) without adding (|157|) . 

This is, in fact, a general phenomenon in the matrix model calculations: for any genus the only 
source for the singular contribution comes from degenerate geometry of curves and is related with 
normalization factor in the matrix integral, that is, the volume of (the orbit of) the unitary group. 
Indeed, the integral itself is a Taylor series in S^s (see formula (4.8) in [2H]), while the unitary group 
volume jnni contributes with the factor 

/Si/n \ 

n n r « =ii G ^/ h ) ( i5§ ) 

i \l=l J i 

where G^ix) is the Barns function, [60 . Now using the asymptotic expansion for the T-function at 
large values of argument and formulas relating the T-function and the Barns function (see, e.g., |61j). 
one finds the asymptotic expansion of the Barns function I(i2j 



N \ Q 2 



S\ „ 1 . 3„ 1„, . ^ B 2h 1 



G 2 (N) = log ( J] r(I) 1 = - logiV - - logiV - -N* + - 2 N\o g 2v + C '(-l) + £ _ N , h _ 2 

(159) 

where B<ih are Bernoulli coefficients and £(s) is the Riemann (^-function. Thus, one obtains that the 
singular contribution in genus h is (25] 

n-ttjjfafi** (160) 

This is the simplest way to pick up the singular contribution, although it can be also done, genus by 
genus, by direct solving the loop equations, like it was demonstrated above for the genus one case. 

Now let us turn to another important issue. In sect. 01 we obtained formula (|129|) . expressing the 
third derivatives of J-q through the quantities 4>i a determined in l|128[) . Since one often interprets T 
as the free energy of a topological string theory, one could naturally associate the third derivative of 
J-q with the tree three-point function in this theory, i.e. represent (|129|l as three "propagators" <pi a 
ending at the same "3- vertex" : 



In such case, one has to associate T\ with the one-loop diagram in this topological theory, i.e. with 
the propagator determinant detf t0l (j)j a : 




Calculating this determinant stems actually to calculating the polynomial determinant det/^ Hj(n a ). 
We already saw that, due to normalization conditions (|57|) and (|74|) . the polynomials Hk(X) corre- 
sponding to the variables tk with k > always have the coefficient k at the highest term \ n ~ l + k while 
the polynomial Hq(X) starts with unit coefficient at A n_1 . Passing from Hi(X), corresponding to the 
variables Si, i = 1, . . . , n — 1, to the basis of monomials A* -1 , one obtains that the total determinant 
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is then (up to a trivial factor n!) just the total Vandermonde determinant divided by the determinant 
of the transition matrix a (|1U4|) . The complete answer for the determinant of 4>j a is then 




(161) 



and the second of matching conditions ()134|) stems now to the condition of nontriviality of det a. 24 

Comparing (|161|) and (j!55|) . one finds that powers of deter and of the Vandermonde determinant 
A(/z) in these expressions perfectly match, i.e. T\ is indeed proportional to the determinant, up to 
non- universal pieces containing M and (arbitrary function of) Si. These pieces remain due to the 
freedom in defining the measure in the path integral. 

Therefore, we conjecture an existence of a diagram technique for calculating the higher genera 
free energy and/or generating function for the correlators 25 . Would such a diagram technique be 
constructed in full, it opens a possibility of calculating the higher genera/multi-point contributions 
in a rather effective way. Therefore, it would be of great practical use to make further checks of the 
conjecture. 



4.4 Relation to topological B-model 

The authors of 23 proposed an anzatz for T\ in the two-cut case (with absent double points). Their 
formula in fact comes from the correspondence between the so called topological B-model on the 
local Calabi-Yau geometry II and the cubic matrix model conjectured in [3]. However, this does not 
completely fix the formula for Ti, leaving room for a certain holomorphic ambiguity, which was fixed 
in |23| basically by some simplicity arguments. 

First of all, introduce the quantities /i7 2 = {^2 — Mi, tH ~~ ^3}; i- e - complexified lengths of the two 
cuts on hyperelliptic plane, and {Si, S2 = S n = to — Si}. Then, one expects 

A(/i) 2 / 3 (/ii + /i 2 -//3-At4) _1 ^ • (162) 

This formula was, indeed, checked for a few first terms of expansion in Si's |23j and it is proven by 
the direct calculation in |64| . 

Below we propose a similar formula for the case of any number of cuts (see also the details in 64 ). 
Let us divide all the branching points into two ordered sets {/x^}™ =1 and {/A }™ =1 and perform then 
a linear orthogonal transformation of fJ-f'^ to the quantities {/^}™ =1 and {^~}™ =1 by 

V$ = l if±$\ (163) 

Taking now n — 1 canonical variables Si, the variable S n = to — J27=i P lower times tk, k = 1, . . . ,p 
(0 < p < n), and choosing an arbitrary set of n + p branching points jj, a . , j = 1, . . . , n + p, following 

24 This is easy to prove. Would be det<r = 0, one obtains that there must exist a polynomial P(X) of degree less or 
equal n — 2 such that 

C 21 ^ = for ; , „ ,. 

This necessarily implies that P(A) has at least one zero at each of the intervals (A2i_i,A2i); otherwise the combination 
under the integral sign is sign definite and the integral cannot vanish. The polynomial P(X) must then have at least 
n — 1 zero and, having the degree not exceeding n — 2, must therefore vanish. 

25 Note that this conjectured diagram technique is different from that of 1321 . Note also that recent paper |63|. where 
the diagram technique of IM2I was extended from calculating resolvents to the free energy calculations, possesses a clear 
disadvantage: intermediate patterns appeared are manifestly non-symmetric w.r.t. field propagators. 
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the same logic as for (|161(l (see also 1)126(1 ). we obtain 



det 



9{Si, S n , t^} 



A(/i a .) • (deter) 



-i 




(164) 



with the same matrix a for any choice of the set of indices {aj}^f and any number p of canonical 

times tk (but only for < p < n). Set all Ma = 1; the Vandermonde determinant A(/j, aj ) then 

combines with the rational factors in the denominator to produce (— 1)"j =1 aj A(/Z^j)/A(/i), where 
A(/U^j) is the Vandermonde determinant for the supplementary set of n — p branching points not 
entering the set {// a ^}"Zf ) whereas A{fi) is the total Vandermonde determinant. In particular, when 
p = 0, splitting /i a as in (|163|) and using formulas (|155() and 1)164)1 . we have 



1 



(i)- 



log det 



d{Si,S n } 



(165) 



where the additional Vandermonde determinant is taken w.r.t. the supplementary variables fij 
the two-cut case it reproduces (f 162)) . 



In 
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